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Ragwitz and Kantz pj propose a correction to a 
method for the reconstruction of Fokker-Planck equa- 
tions from time series data. In §, |, |, § a method 
was presented which directly applied the mathematical 
definitions of the drift and diffusion term || 
for an estimate from time series. Here different moments 
of conditional probability densities (pdf) for finite step 
sizes A in the limit A — > have to be estimated. Rag- 
witz and Kantz state that previous results have not been 
checked and that indispensable finite time step A correc- 
tion have to be employed for reliable estimates of D^ 2 > . 
We want to add the following comments. 

Ragwitz and Kantz base their investigation on an es- 
timate of the finite time conditional probability in terms 
of a Gaussian, eq. (7) of their paper. There is, however, 
no reason that for finite A the conditional pdf is Gaus- 
sian. The exact expressions for the conditional moments 
up to the order A 2 can be unambigously derived from 
the Fokker-Planck equation |Q : 

(x-x \x Q ) = A£>« + (1) 

+ -A 2 (£>«)' + D^(D^)"] + 0(A 3 ) 

{(x - zoflxo) = 2AD^ + A 2 [(£«) 2 + (2) 
+ 2D^(D^)' + DW(DW)' + dW(dW)"} + 0(A 3 ), 

For ((x — xq) 2 \xo) the ansatz (7) of [Q neglects the last 
two terms which are important for processes involving 
multiplicative noise. However, intermittency of turbu- 
lence is related to a multiplicatice process. This remark 
especially applies to the wind data presented in |l| . The 
validation of their method based on a Langevin (equation 
(9)) only works since purely additive noise is considered. 

For the investigation of turbulence Ragwitz and Kantz 
claim to obtain remarkable correction, as shown in their 
Fig. 6. Our approach in H is based on an estimate of the 
diffusion term using the limit A — > yielding a depen- 
dency which can be approximated by a low order polyno- 
mial. In order to improve on this estimate the coefficients 
of this polynomial have been varied in such a way that 
the solution of the corresponding Fokker-Plank equation 
for finite values of A yields an accurate representation 
of the measured one. In other words, in a second step, 
we have performed a parametric estimation of the diffu- 
sion term. In Figure 1 we present a case where a large 
correction to the A — > estimation of had to be 
introduced (usually corrections are of the order of some 




A//, 



FIG. 1: A dependence of D {2) ( u = a 00 ,r = L/2, A) /cr^ for 
different correction terms; squares without correction term, 
circles with ((D' 1 ') 2 ), and crosses with the correction term of 
L denotes the integral length, (T^ the rms of the velocity 
increments at large scales. Only for A/l mar > 1 Markovian 
properties hold, and estimations of are senseful. The 

optimal value of D' 2 ' based on verifications is indicated by 
an arrow, further detail see H. 

percentages). For finite values of A the estimated values 
of D^ 2 ' clearly differ from the limiting case A — > 0. Es- 
timations with different "correction" terms for and 
for finite A values may fake large corrections values. Tak- 
ing the limit A — > these deviations vanish within the 
error. 

The range of A which can be taken for the estimate 
of D^ 1 must be choosen carefully in order to ensure that 
the Markovian property holds, see f| . Since for each es- 
timated value of Z}( 2 ) a finite number of data point is 
used, a statistical error can be estimated for (cf. 
Fig. 1). These errors naturally increase considerably for 
large values of x (compare Fig. 6 [Q and Fig 13 M.) To 
conclude, a deeper understanding of finite time correla- 
tions are of interest. The correct terms of higher orders 
in A may be used to improve on the estimation of 
drift and diffusion terms. Up to now the best way for 
estimating diffusion coefficients is to combine a nonpara- 
metric estimate for A — > with a functional ansatz, i.e. 
a suitable polynomial ansatz. Refining the estimates of 
the coefficients by parametric methods as for instance de- 
scribed by Q leads to improved results by a comparison 
of measured and calculated conditional pdfs at finite A. 
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